Hadron yields from thermalized minijets 
at RHIC and LHC 
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Abstract 

We calculate the yields of pions, kaons, and (/)-mesons for RHIC and LHC energies 
assuming thcrmodynamical equilibration of the produced minijets, and using as 
input results from pQCD for the energy densities at midrapidity. In the calculation 
of the production of partons and of transverse energy one has to account for nuclear 
shadowing. By using two parametrizations for the gluon shadowing one derives 
energy densities differing strongly in magnitude. In this publication we link those 
perturbatively calculated energy densities of partons via entropy conservation in an 
ideal fluid to the hadron multiplicities at chemical freeze-out. 
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1. Introduction 

Particle production in high-energy heavy-ion 
reactions at the BNL-RHIC and CERN-LHC 
colliders will soon provide interesting insight into 
nuclear modifications of semi-hard processes [Q. 
This is because pQCD processes involving gluons 
in the initial state may dominate the inelastic AA 
cross-section at collider energies. In particular, they 
might lead to a better understanding of the gluon 
distribution in large nuclei, which is not accessible 
in DIS. 

In Q the effect of nuclear shadowing of 
the parton distribution functions on the charged 
particle multiplicity at midrapidity has been 
investigated assuming no rescattering between the 
produced minijets and the hadrons they fragment 
into. Here, we will take the opposite point of view 
and assume maximal rescattering, i.e. local thermal 
and chemical equilibrium of the minijets. We 
compute final-state hadron multiplicities of various 
hadron species under the assumption of entropy 
conservation. 

In ||] we calculated the initial conditions at 
RHIC and LHC by means of pQCD above the 
semihard scale pr — 2 GeV to derive the number 
and energy densities of partons at midrapidity. 
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In that calculation we explicitely included the 
shadowing effect on the parton distribution 
functions entering the formulas for the production 
of flavor / = g,q,q in the minijet approach. We 
employed two different parametrizations for the 
shadowing effect accounting for weak and strong 
gluon shadowing, respectively, |^,|4|. A direct 
consequence of the shadowing effect is the decrease 
in the production of partons of given momentum 
Pt, i.e. a decrease of transverse energy production 
at midrapidity. We calculate the first Et moment 
with and without shadowed pdf's and with a 
cut-off function e{y) ensuring that we only count 
scatterings into the central rapidity region (\y\ < 
0.5): 
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Motivated by the factorization in QCD one can 
assume that the production of transverse energy in 
AA collisions can be split up into a hard and a soft 
contribution as 



\ppi 



(2) 



With an energy independent value of crf^ = 32 
mb one derives g aPP {Erf/ = 15 mb GeV. 
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With TauAu — 29/mb one can derive the soft 
contribution (i.e. the one for < 2 GeV) to the 
energy density for RHIC as Ssoft =33.7 GeV/fm . 
When comparing to the soft contributions one finds 

with ahardiET)^, = i^'+<^''+CT'^)hard{ET)^ard 

that the ratio of soft to hard contribution 

^"'^ ^ — 

is Rsh = 0.47 for no, Rsh = 0.77 for strong, and 
Rsh = 0.47 for weak shadowing, respectively, since 
£ = 71.5 GeV/fm^ for no, e = 43.6 GeV/fm^ for 
strong, and e = 71.9 GeV/fm'^ for weak gluon 
shadowing . This imphes that at RHIC the soft 
component could even dominate if it were energy 
independent and also unaffected by the shadowing 
effect. 

For LHG we only calculated the contribution 
of the gluons that strongly dominate all partonic 
processes due to the large distribution function at 
small momentum fractions. The results for the 
energy densities for no, strong, and weak gluon 
shadowing, then are £g = 1229.7 GeV/fm'^, Eg = 
144.8 GeV/fm^ and = 678.6 GeV/fm^ With 
TpbPb = 32/mb one can derive the energy density 
from the soft part and has Esoft — 35.8 GeV/fm'^ 
which is slightly larger than at RHIC due to 
the larger nuclear overlap function, i.e. the larger 
number of effective scatterings in the Glauber 
picture leading to the transverse energy production. 
The relative weight 

p _ C^soft {Et)s f.s 
(^hard \t^T)h 

between soft and hard contributions therefore 
changes and becomes Rsh — 0.029 for no, Rsh = 
0.25 for strong, and Rsh = 0.052 for weak 
shadowing, respectively. 

Therefore the soft contributions gain much more 
weight in this naive picture due to the strong 
effect of shadowing on the small-x gluons (for 
further details of the calculation, the shadowing 
parametrizations, etc. see Q). 

2. Hadron Multiplicities at Chemical Freeze- 
Out 

Having calculated the energy densities at midrapid- 
ity for RHIC and LHC in pQCD we connect £i to the 
number of hadrons at midrapidity by assuming an 
ideal fluid that is characterized by entropy conserva- 
tion from the quark-gluon plasma to the hadron gas, 
dSi/dy = dSf /dy [||. We relate the energy density 
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Table 1. Hadron yields at freeze- out with initial condi- 
tions from pQCD for LHC. The initial temperatures are 
calculated from the energy densities produced via hard 
processes and are for a three flavor quark-gluon plasma 
With niu ~ rud = MeV and rUs = 150 MeV, respec- 
tively. 

of the quark-gluon plasma to its entropy density via 
the bag model equation of state 0] . We account for 
u, d, s quarks (with masses ruu = md = 0, nis = 150 
MeV), the antiquarks, and gluons. The total pro- 
duced entropy dSi/dy is obtained from the entropy 
density at midrapidity as dSi/dy = ViSi with the 
initial volume of the central region Vi ~ 7ri?^r,;, 
which is numerically Vi = 12.9 fm^ for Au -\- Au 
and Vi = 13.4 fm^ for Pb + Pb at b = with 
r = 0.1 fm/c. Since we assumed an ideal fluid the 
total entropy is conserved throughout the expan- 
sion until freeze-out which is chosen here to hap- 
pen at a temperature Tpo = 160 MeV. For sim- 
plicity we furthermore assume vanishing chemical 
potentials in the central rapidity region, i.e. that 
all conserved currents are identically zero. If this 
were not true one would have to multiply by factors 
exp{^i/T) (in Boltzmann approximation). The en- 
tropy density of the hadronic fluid is calculated as- 
suming an ideal gas composed of all hadrons up to 
a rest-mass of 2 GeV. Their respective occupation 
numbers are given by Fermi-Dirac or Bose-Einstein 
distribution functions, respectively. Thus, Tpo and 
dSf / dy = dSi/dy determine the multiplicity of each 
hadron species uniquely |8[^. Feeding from post 
freeze-out decays of heavier resonances is also taken 
into account. 

3. Results 

With the model outlined above and the energy 
densities derived above we calculated the number of 
a variety of hadrons at midrapidity. We also include 
the multiplicities due to the soft contributions and 
quote the initial temperatures for a QGP of three 
flavors. For LHC we derived the yields shown in 
table |l| and for RHIC the ones in table 0. For 
the latter one, one can clearly see that there is 
no change in the hadron yield for weak shadowing 
and the unshadowed case, respectively, due to the 
almost identical energy density serving as an input 
for the calculation. 
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Table 2. Hadron yields at freeze-out with initial 
conditions from pQCD for RHIC. 



4. Conclusions and outlook 

We computed the rapidity densities of a variety 
of hadrons based on the assumption of entropy 
conservation of an ideal fluid. The entropy densities 
were derived from the energy densities which in turn 
were caluclated by means of pQCD |^ . We find that 
in the hmit that the minijets equihbrate locaUy the 
effect of shadowing on the hadron yield is not as 
large as on the pure partonic degrees of freedom. 
This can be seen e.g. in the ratio of energy densities 
between unshadowed and strongly shadowed gluons 
at LHC, which is about a factor of 9, while the ratio 
of the hadron yields is only a factor of 5. Since 
vanishing net baryon and strangeness densities were 
assumed, the relative depletion of shadowed to 
unshadowed gluon distribution is independent of 
the particle species. For RHIC, even without 
shadowing we only get about 300 pions since the 
perturbative calculation, entering via the energy 
densities, was performed with the cut-off pt — '2. 
GeV. Therefore, the soft contribution constitutes 
a significant part of the transverse enerry 
and therefore of the particle multiplicities 
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